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Abstract 

We present synthetic spectral fits of the typical Type lb SN 1999dn and 
the Hydrogen Rich lb SN 2000H using the generalized non-local thermodynamic 
equilibrium stellar atmospheres code PHOENIX. We fit model spectra to five epochs 
of SN 1999dn ranging from ten days pre-maximum light to 17 days post-maximum 
light and the two earliest epochs of SN 2000H available, maximum light and six 
days post-maximum. Our goal is to investigate the possibility of hydrogen in 
Type lb Supernovae (SNe lb), specifically a feature around 6200A which has 
previously been attributed to high velocity Ha. In earlier work on SN 1999dn we 
found the most plausible alternative to Ha to be a blend of Si II and Fe II lines 
which can be adjusted to fit by increasing the metallicity. Our models are simple; 
they assume a powerlaw density profile with radius, homologous expansion, and 
solar compositions. The helium core is produced "burning" 4H^He in order to 
conserve nucleon number. For models with hydrogen the outer skin of the model 
consists of a shell of solar composition. The hydrogen mass of the standard 
solar composition shell is Mh ^ IO'^Mq in SN 1999dn and Mh ^ 0.2Mq for 
SN 2000H. Our models fit the observed spectra reasonably well, successfully 
reproducing most features including the characteristic He I absorptions. The 
hydrogen feature in SN 1999dn is clear, but much more pronounced in SN 2000H. 
We discuss a possible evolutionary scenario that accounts for the dichotomy in 
the hydrogen shell mass between these two supernovae. 

Subject headings: supernovae: 1999dn, 2000H 

1. Introduction 

1.1. Supernovae Nomenclature 

Supernovae (SNe) are divided into two main subtypes. Type I and Type II, based on 
the presence of hydrogen in their spectra. SNe of Type I are then further subdivided into 
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Types la, lb, and Ic (IWheeler &: Harknesslll990l ). SNe la are characterized by strong Si II 
absorption near 6IOOA and a characteristic sulfur 'W shaped feature at 5400A. SNe lb are 
classified by the presence of optical He I lines in their spectra, and SNe Ic are noted for their 
lack of helium and silicon. SNe la are thought to be caused by the thermonuclear explosion 
of a Chandrasekhar mass white dwarf star while SNe II, lb, and Ic are caused by the core 
collapse of a massive star. The general consensus is that most SNe Ib/c result from mass 
loss due to interaction with a binary companion and that SNe l ib also result from binary 
interaction (jPodsiadlowski et al.lll992l . ll993uNomoto et al.lll995uHeger et al.ll2003l ). 

Supernovae are classified by their spectra. As early spectra become more commonly 
available, trends in the variation of certain SNe types over the course of time have emerged. 
This is apparent in SN 1993 J which was originally classified as a SN II for its conspicuous 
hydrogen but later showed helium absorptions characteristic of a SN lb. SN 1993J is now 
classified as a SN lib which is believed to be the explosion of a s tar that has undergori e 
enough ma ss loss to lose a l arge portion of its hydrogen envelope (Filippenko et al.lll993l ). 
SN 1987K flFilippenkolll988h and SN 2008ax dchornock et al.lboioh are other examples of 
this transitory class of SNe. This explains the distinct Balmer lines seen at early times which 
give way to significant helium feature s as tim e progresses. This type of transition-like object 
was first suggested by lWoosley et al.l ( 119871 ). It has been suggested that if mass loss from a 
SN II progenitor can produce SNe lib then an even larger mass loss woul d result in furthe r 
stripping of the hydrogen envelope and would be characteristic of SNe lb (lFilippenkdll988l ). 
This suggests that perhaps we may see hydrogen in ea rly spectra of SNe lb if there was a 
large enough mass of hydrogen left prior to explosion. iNomoto. Iwamoto. &: Suzuki! (jl995[ ) 
also suggest several binary mass loss scenarios that might lead to a spectroscopic sequence of 
SN types from Iln— )-IIb— >Ib— )-Ic. By understanding the behavior of these transitory types of 
SNe we are able to better extrapolate back to the state of the progenitor star. We therefore 
investigate the early and late time spectra of a typical Type lb S N 1999dn and the early 
spectra of a Hydrogen Rich SNe lb SN 2000H (IBranch et al.l 120021 ) in search of clues as to 
whether or not a thin hydrogen skin exists prior to the supernova explosion. 



1.2. Prior Investigations 



The idea that there exists an underlying physical connection between SNe II and SNe lb 
is not a new one. As more SNe are discovered at earlier times we see a number of these 
SNe lib, wh ich appear to be trans ition-like objects from Type Il^^Ib. The recent discovery of 
SN 2008ax (jChornock et al.ll2010l ) h as led to more de t ailed investigations into the connection 
between SNe lib and SNe lb, and IChornock et al.l (120101 ) suggested that a full non-local 
thermodynamic equilibrium (NLTE) treatment is necessary for a conclusive analysis of the 
details of the hydrogen skin. 
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Branch et al. 



20021) conducted a direct analysis of a selection of SN lb spectra obtained 



by lMatheson et al.l ( 1200 ih using the supernova spectrum synthesis code SYNOW. They con 
eluded that a small amount of hydrogen is present in perhaps all SNe lb spectra, specifically 
SN 1999dn, and that even m ore is seen in what they term Hydrogen Rich SNe lb: 2000H , 



1999di, and 1954A. Others flElmhamdi et al.l 120061 : lOeng et al.l l2000l : IParrent et al.l l2007f ) 



have also used SYNOW to investigate the presence of hydrogen in SNe lb. 

Because the relative simplicity of SYNOW a first principles non-local thermodynamic 
equilibrium (NLTE) treatment is needed to analyze the presence of hydrogen in detail. While 
a SYNOW analysis can determine whether a particular ion is there or not and can also 
indicate the relative strength and velocity extent of that species, one cannot use SYNOW 
to determine quantative abundan ces. We use the NLTE stellar atmospheres code PHOENIX 
(IHauschildt Sz BaronI |2004J . Il999l ) to accomplish this. We perform a time-series analysis of 
the spectra of SN 1999dn, a typical SN lb, as wel l as some calculations of early time spectra 
of another SN lb, SN 2000H. iBenetti et all koO(h originally classified SN 2000H as a SN lib 
due to strong and possible if^ absorptions. 

We wish to accomplish the following: (1) Produce hydrogen and hydrogen-free fits to 
relevant epochs of SN 1999dn and SN 2000H; (2) Constrain the possible ejection velocity 
and mass of the hydrogen skin. 



2. Methods 
2.1. General Methods 

We generate synthetic spectra of SNe lb using the non-local thermodyna mic equilibrium 



(NLTE) stellar atmospheres code PHOENIX (jHauschildt &: BaronI l2004i Il999l . and references 



therein) The inclusion of NLTE calculations allows us to investigate the importance of non 
thermal processes in reproducing SN lb spectra. We explore the possibility of the existence 



of hy drogen in SNe lb by adding a thin skin of solar composition material (iGrevesse et al. 



20071 ) to a helium core (which also has solar compositions of metals). The helium core is 
constructed by beginning with material at solar composition and then "burning" 4H— >-He. 
We then add the hydrogen skin with solar abundances by replacing the abundances above a 
specified velocity, Vh, which is varied to find the best fit without disturbing the helium core 
interior to Vh- Our models extend out to a maximum velocity of 21000 km s~^ regardless 
of whether they are pure helium cores or include a hydrogen skin. We approximate the 
gamma-ray deposition function by a constant parameter that is adjusted for each epoch. 
The assumption of a sharp boundary between the hydrogen and helium layers is not realistic, 
but is simple and eliminates the need to add extra parameters. This approximation should 
not significantly effect our conclusions. 
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The earliest model of SN 1999dn is only M = O.SSM© as compared with all the other 
days that are M = 38Mq. This difference in mass was necessary in order deal with the large 
variation in the dynamic range of the atmosphere. Nevertheless our models are internally 
consistent. The very large mass of the progenitor is due to the straightforward extrapolation 
of the steep powerlaw density profile to low radii. Here it is just a calculational expedient 
and does not imply that the progenitor was anywhere close to this massive. 



3. Results 
3.1. SN 1999dn 



Fig, m shows the observed spectra for several epochs of SN 1999dn ( iMatheson et al. 



2001uDeng et al.ll2000l ). While we seek the overall best fit to the spectra, emphasis is placed 
on the absorption feature around 6200 A which fades after maximum light. The observed 
spectra are smoothed with a 10 point boxcar average in order to remove some noise. 

We explored variations in multiple parameters for each epoch adjusting the total bolo- 
metric luminosity in the observer's frame (specified by a temperature, Tmodei)) the photo- 
spheric velocity, the metallicity, the density profile, gamma-ray deposition, and hydrogen 
mass. By varying these parameters we have arrived at a set of parameters which we con- 
sider the best and we follow that model through each epoch. The relevant parameters for 
SN 1999dn are given in Table □ and in Table [2] for SN 2000H. The results of the NLTE 
PHOENIX calculations are shown in Figs. |2H8l We treat the following species in NLTE: H I, 
He I-II, C I-III, N I-III, O I-III, Ne I, Na I-II, Mg I-III, Si I-III, Ca I-III, and Fe I-III. In all 
the spectra we successfully reproduce the characteristic He I absorptions at 5876A, 6678A, 
and 7065A. 



Fig. [2] shows the the early epoch Aug. 21 for SN 1999dn flMatheson et all l200lL -10 
days) the synthetic spectrum with hydrogen at 19000 km s^^ fits the 6200A feature better 
than a pure helium core. In the maximum light spectrum, Aug. 31, the hydrogen feature 
at 19000 km s^^ does not properly blend with the feature just to the red (Fig. [3]), but this 
is likely due to our assumption of a sharp abundance boundary. The detached-like profile 
(blueshifted absorption with fiat topped emission) is a consequence of the abrupt cutoff of 
our model. The hydrogen velocity is 19000 km s~^ and the outer velocity layer of our model 
is 21000 km s~^. A smoother variation in the abundances would reduce the sharp cutoff that 
leads to the detached profile. 

By Sept. 10 the hydrogen feature has faded in both the observed spectrum as well as 
the synthetic spectra as seen in Fig. HI This is reinforced by the lack of a feature in the later 
days of the observed spectra. 
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In Fig. [5] and Fig. [6l on Sept. 14 and Sept. 17 respectively, we find that the inclusion 
of hydrogen has little if any effect on the spectrum and so conclude that the thin skin of 
hydrogen is no longer visible due to geometric dilution. 

The He I absorptions in many of the spectra appear to be slightly to the red of the 
observed absorptions. This is probably due to the assumption of a constant gamma-ray 
deposition parameter. By increasing the amount of gamma-ray deposition we increase the 
non-thermal components of the spectrum which increases the depth of the He I absorptions 
and also the speed of highly non-thermal lines. A more thorough treatment of the nickel 
distribution would likely improve this, but examining the exact amount of nickel mixing is 
beyond the scope of the present work. 



3.2. SN 2000H 



The earliest spectrum available of SN 2000H is at maximum light, Apr 21 2000 (IBarbon et al 



20091 ). Fig. [7] shows a very deep and wide absorption trough in the region of absorption. 
Part of the absorption feature could be due to blending with other elements due to circum- 
stellar interaction. Towards the end of this section we discuss the narrow Na D absorption 
near 5800A and the information we can gather about the circumstellar medium (CSM) from 
it. 

Part of the line profile could be due to a low velocity CSM region but the depth of the 
line would suggest that we should also see evidence of Hp, which is lack ing. This supernova 
is highly reddened. We adopt a reddening value of E{B — V^) = 0.61 (IBenetti et al.l 120001 ) 
and a recession velocity of Vrec = 3945 km s~^ taken from the NEE0 database. 



We find that the spectrum of SN 2000H at +6 days post max light (IBarbon et al.ll2009l ) 
can also be well reproduced by adding hydrogen to our models as is shown in Fig. EJ We 
find that a hydrogen velocity of vh = 15000 km s~^ is necessary to fit the hydrogen feature 
around 6200A. The amount of hydrogen included in the 2000H model {Mh <^ 0.2Mq) is 
significantly greater than in our SN 1999dn models , (Mh <; lO'^M^), which is consistent 
with the classification of a Hydrogen Rich SN lb by [Branch et al.l (|2002[ ). This is necessary 
to fit the hydrogen feature that persists to as late as +6 days post-max light. A larger 
amount of hydrogen suggests that SN 2000H is c loser to the Type lib subtype due to a lower 
degree of envelope str ipping than in SN 1999dn (IBenetti et al.ll2000l : iBlondin fc Tonryll2007l : 



Chornock et al. 



20101 ). 



Figs. lUHTU] show two of the observed spectra of SN 2000H both with and without a 20 
point boxcar smoothing. The narrow Na D absorption around 5800A appears to fade quickly 
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since it is not seen in the spectrum taken 6 days later. In Fig[TT]we present a comparison of 
a SYNOW spectrum of only Na I constrained to 3000-5000 km s~^ as well as the case where 
the Na I is detached from the photosphere but constrained to 3000-3500 km s~^ with the 
maximum light spectrum of SN 2000H with 20, 10, and 5 point boxcar smoothing. The Na 
I and the 20 point boxcar spectrum look like the best fit, but since the SYNOW spectrum 
appears too broad to fit the 5 point boxcar the detached spectrum is a better fit in that case. 
In either case one finds a characteristic velocity of the Na D line to be about 3000 km s~^. 

We believe the Na D feature to be formed fro m the surroundin g GSM. The SYNOW 



velocity is large compared to typical wind velocities. lYoon et al.l (120101 ) found a wind velocity 
for different SN Ib/c progenitors to range from ~ 200 km s~^ to ~ 2400 km s~^ which at 
the upper limit could explain the high velocity we find for the Na D feature from the GSM, 
however with a hydrogen shell of ~ 0.1 Mq, the wind is more likely to be of order 500 km s^^. 
The high velocity could also be due to radiative acceleration at shock breakout. We find that 
the total energy required to accelerate this mass to be i^^in = 2 x 10^^ ergs(p/10~^^ g cm^'^) 
where the density is obtained from the SYNOW fit assuming an ionization fraction of 10^^. 



4. Discussion 

We find that the results from the PHOENIX calculations suggest a hydrogen mass of 
Mh <i, 10~^M ( 7i for S N 1999dn. This corresponds to a hydrogen velocity vh ^ 19000 km s~^. 



Tanaka et al.l ( 12009I ) placed a limit of ^ 5 x W^Mq on the hydrogen mass in SN 2008D 



and a limit on the Doppler velocity of approximately vh ~ 18500 km s ^. Others have also 



suggested that hydrogen is responsible for this feature around 6200A. IPeng et al.l (2000|) 
suggested a hydrogen velocity for SN 1999dn of 19000 km s^^ on Aug. 21, and 18000 km s~^ 
on A ug. 31, and tha t the feature becomes blended and then taken over by G II in later 



days. iKetchum et al.l (120081 ) investigated G II, Si II, Fe II, and Ne I as possible candidates 
as well but found that without increasing the metallicity and abundances beyond physically 
reasonable amounts the 6200A feature was not easily reproduced in early epochs. This 
seems to be consistent with our calculations, reinforcing the the existence of hydrogen in 
early spectra of SNe lb. 

We also find that the results from the PHOENIX calculations suggest a hydrogen mass 
of Mh <s O.2M0 for SN 2000II. This corresponds to a hydrogen velocity in our models of 
Vh ~ 15000 km s~^. This mass estimate is likely an upper limit since we assume a single 
power law density profile for both the helium core and hydrogen shell. The width of the 
6200A feature early on is likely due to some blending with some circumstellar material since 
it disappears in the spectrum 6 days later and the hydrogen feature is more clear. We know 
that circumstellar interaction was taking place at maximum light from the presence of the 
narrow Na D feature that quickly fades. The larger amount of hydrogen is more characteristic 
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of that seen in a SN lib such as SN 1993J or SN 2008ax with mass estimates on the order 



of Mh ^ O.IMq flWheeler et al.lll994J : IChornock et al.lbOlOf ). 



The possibihty of a hyd rogen skin existing prior to the supernova explosion is interesting. 
Podsiadlowski et aL (Il992) first discussed the evolution of binary systems that would lead 
to mass stripping. iNomoto et al.l ( 1l995h discuss a continuum of supernovae types from 
Iln— >^IIb— J-Ib^lc which would be consistent with the idea of a varying degree of hydrogen 
mass loss due to a strong wind or binary interaction in a non-conservative mass transfer 
scenario. The energy rel e ased i n a binary merger would be sufficient to eject some of the 
envelope. I Wheeler et al.l (jl994j ) suggested that if hydrogen is present in SNe lb then the 
amount of hydrogen in SNe lb is less than that found in the Type lib SN 1993J, which was 
about 0.1 — 0.5 Mq (IFilippenko et al.lll993l ). This is consistent with ou r results of SN 1999dn 
but su ggests that SN 2000H is more accurately classified as a SN lib. I Chevalier Sz Fransson 
(|2006[ ) also suggest that there is a continuous distribution of SNe between SNe lib where the 
hydrogen is clearly present and SNe lb where the hydrogen is quite weak. This continuous 
distribution is likely to be the explanation for the sequence of SNe between type lib and lb 
and the varied amounts of envelope stripping are a result of binary interaction and/or strong 
stellar winds. SN 1993J is known to have been in a binary and as a r esult had some of its 



hydrogen envelope lost du e to the merging of the two stars in the binary (jPodsiadlowski et al. 



1993 



Nomoto et al.lll995h 



Recently lYoon et al.l (l2010[ l studied binary progenitors of SNe Ib/c including the effects 
of rotation and various reduced wind-loss formulations. They find that even in the case of 
Case A and B mass transfer that the mass transfer is non-conservative and that a limited 
range of helium cores (3.5 ^ M <^ 4.5 Mq, for solar metallicity and 3.5 ^ M 8 Mq, for 
Zq/20) that the wind loss in the helium WR stage is low enough that a hydrogen skin of 
the order of Mh ~ 10^^ — 10^'^ can be retained. 

Even though there is a rather continuous distinction between SNe lib and hydrogen 
rich SNe lb, the total hydrogen mass in the envelope is rather similar Mh ~ 0.1 — 0.3 Mq. 
These objects are likely from binary progenitors with stable Case C mass transfer, which 
leads the donor to shrink below the Ro c he limit wh en the envelope reaches a critical mass 
of about 0.3 Mq ( jPodsiadlowski et al.l Il992l Il993l ). Further reduction of the envelope 
mass then follows via ablation o f the hydrogen envelope from the wind of the secondary 
( jPodsiadlowski et al.l Il992l . Il993l ). However, studies of SNe lb suggest that many if not 
most have a much thinner skin of hydro gen similar to what we find here for SN 1999dn 
( jBranch et al.ll2002[ IChornock et al.ll2010l ). This could likely be explai ned by a binary sce- 



nario that leads to a co mmon envelope that is ejecte d from the system (jPodsiadlowski et al. 



19921 ). As described in iPodsiadlowski et al.l (119921 ) this will lead to either a cataclysmic 
variable (if the core mass is less than 1.4 Mq) or to a low mass helium core, which would 
explode as a SN lb, but in the common envelope environment it would not be hard to keep 
a thin hydrogen skin. If this is a common channel for SNe lb formation than one would 
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expect many to have thin hydrogen envelopes which will be visible in early spectra. Since 
massive helium cores will blow strong winds that would quickly remove a 10^"^ Mq envelope 
this lends support that many SNe lb progenitors are the result of relatively low mass helium 
cores. 



5. Conclusion 

Our results indicate that hydrogen does exist in the Type lb SN 1999dn and SN 2000H 
with an upper mass limit of Mh ^ 10~^Mq and Mh ^ 10~^Mq respectively. The existence of 
such a small amount of hydrogen is to be expected as suggested by the continuous spectrum 
of SN types from II to Ic. Further analysis of a more robust sample of SN lb is needed for a 
definitive claim, but the existence of a hydrogen skin in SNe lb is seen in the observed spectra 
as well as supported by synthetic calculations. This implies that the typing of core-collapse 
SNe is more useful as a means of determining the state of the ejecta at a given epoch as 
opposed to inferring the state of the progenitor at the time of explosion. While it is still 
useful to determine a SN type we need to set a uniform epoch, such as maximum light or 
later, in order to be consistent and conclusive when it comes to classification. 

We thank Phillip Podsiadlowski for extensive discussions and tutelage on the nature 
of binary interactions in stripped envelope supernovae. This work was supported in part 
NSF grant AST-0707704, and US DOE Grant DE-FG02-07ER41517. This research used 
resources of the National Energy Research Scientific Computing Center (NERSC), which is 
supported by the Office of Science of the U.S. Department of Energy under Contract No. 
DE-AC02-05CH11231. 
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Table 1. Model Parameters: 1999dn 



Epoch 


2^model (K) 


vo (km s ^) 


vh (km s ^) 


Mh (Mo) 


-10 


6000 


11000 


19000 


2.95 X 10-4 





5250 


10000 


19000 


2.94 X 10-3 


10 


5400 


7000 


19000 


2.63 X 10-4 


14 


5000 


7000 


19000 


2.63 X 10-4 


17 


4600 


7000 


19000 


3.38 X 10-4 



Note. — Relevant Model Parameters for SN 1999dn 



Table 2. Model Parameters: 2000H 



Epoch 


Tmodel (K) 


vq (km s ^) 


Vh (km s ^) 


Mh (Mq) 





6000 


11000 


15000 


0.124 


+6 


5600 


10000 


15000 


0.125 



Note. — 



Relevant Model Parameters for SN 2000H 
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^ 1.00 

CO 

'c 
CO 

^ 0.10 



0.01 




I I I . 



3000 4000 5000 6000 7000 8000 900010000 

A (A) 



Fig. 1. — SN 1999dn Epochs: Emphasis i s placed on the 6200A absorption featur e that fades 
after maximum hght. Observed spectra (IMatheson et aLll2001uDeng et aLll2000[ ) have been 
smoothed with a 10 point boxcar. These spectra were all obtained from the Supernovae 
Spectrum Repository (SUSPECT) http: / /suspect. nhn.ou.edUj maintained by the University 
of Oklahoma. 
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Fig. 2. — SN 1999dn (-10) Observed spectrum ( jPeng et al.ll2000l ) smoothed with a 10 point 
boxcar compared to synthetic spectra with and without hydrogen at Vh = 19000 km s~^. 
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Fig. 3. — SN 1999dn (0) Observed spectrum (jPeng et al.l l2000l ) smoothed with a 10 point 
boxcar compared to synthetic spectra with and without hydrogen at Vh = 19000 km s~^. 
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Fig. 4. — SN 1999dn (+10) Observed spectrum (IMatheson et al.l l200lh smoothed with 
a 10 point boxcar compared to synthetic spectra with and without hydrogen at vh = 
19000 km s-\ 
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Fig. 5. — SN 1999dn (+14) Observed spectrum ( jPeng et al.l 120001 ) smoothed with a 8 point 
boxcar compared to synthetic spectra with and without hydrogen at Vh = 19000 km s~^. 



- 16 - 



(+17) observed 
Hydrogen free 
Hydrogen @ 1 9000(km/s) 




4000 6000 8000 10000 
A (A) 



Fig. 6. — SN 1999dn (+17) Observed spectrum (IMatheson et al.l 1200 ll ) smoothed with 
a 10 point boxcar compared to synthetic spectra with and without hydrogen at Vh = 
19000 km s-^ 
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Fig. 7.— SN 2000H (0) Observed spectrum ( IBarbon et all 120091 ) smoothed with a 20 point 
boxcar compared to synthetic spectra with and without hydrogen at Vh = 15000 km s~^. 
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Fig. 8. — SN 2000H (+6) Observed spectrum ( iBarbon et al.ll2009[ ) smoothed with a 20 point 
boxcar compared to synthetic spectra with and without hydrogen at Vh = 15000 km s~^. 
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Fig. 10.— SN 2000H epochs flBarbon et al, 



20091 ) smoothed with a 20 point boxcar. 
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Fig. 11.— SN 2000H (0) f lBarbon et al.ll2nn9h with a 20, 10, and 5 point boxcar plotted with 
a SYNOW spectrum of Na I between 3000-5000 km and detached from the photosphere 
but constrained to a velocity between 3000-3500 km s~^. 



